Localized magnetic states in biased bilayer and trilayer graphene 
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Abstract 

We study the localized magnetic states of impurity in biased bilayer and trilayer graphene. It is 
found that the magnetic boundary for bilayer and trilayer graphene presents the mixing features of 
Dirac and conventional fermion. For zero gate bias, as the impurity energy approaches the Dirac 
point, the impurity magnetization region diminishes for bilayer and trilayer graphene. When a 
gate bias is applied, the dependence of impurity magnetic states on the impurity energy exhibits 
a different behavior for bilayer and trilayer graphene due to the opening of a gap between the 
valence and the conduction band in the bilayer graphene with the gate bias applied. The magnetic 
moment and the corresponding magnetic transition of the impurity in bilayer graphene are also 
investigated. 

PACS numbers: 73.20.Hb, 81.05.Uw, 73.21. Ac 
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I. INTRODUCTION 



The intense research currently devoted to graphene, a two-dimensional carbon honeycomb 
lattice, has uncovered a wealth of fascinating properties such as the anomalous quantized 
Hall effect, the absence of the weak localization and existence of the minimal conductivity 
E, [3, 14, 0]. Graphene has a high mobility, its carrier density is controllable by an applied 
gate voltage [2|] and a spin-orbit interaction 

Graphene structures have been the focus of much interest [12 



26|. In particular, adatoms may be positioned on graphene by current nanotechnology[2l|, 




rendering the study and manipulation of local electronic properties. Ab initio calculations 



for transition metal adatoms 
ments. Recently Uchoa et al 



22l| s how a tendency to the formation of local magnetic mo- 
] examined the condition for the emergence of localized 



magnetic moment on adatoms with inner shell electrons on a single layer graphene. It is 
found that the impurity magnetization boundary exhibits anomalous characteristics. In 
contrast to the case of an impurity in an ordinary metal, the impurity can magnetize for 
any small charging energy due to the low density of state(DOS) at the Dirac point. On 
the other hand, detailed experimental studies |26| on multi-layer graphene showed a marked 
modification of the electronic structure with the number of layers. Hence, we expect ll(] a 
qualitative difference in the magnetic properties of the adatoms on multilayer graphene; an 
issue which we address here by inspecting the localized magnetic state of an impurity in a 
biased bilayer and trilayer graphene. We find that the size of the magnetic region decreases 
rapidly compared with that in monolayer graphene, the impurity can magnetize even when 
the energy of the doubly occupied state is below the Fermi level, and the impurity magne- 
tization region is asymmetric due to the special nature of the quasiparticles having mixed 
features of Dirac and conventional fermions. When a gate bias is applied, the dependence 
of the impurity magnetic states on the impurity energy for a bilayer graphene exhibits a 
different behavior from that for a trilayer graphene due to the opening of a gate-induced 
gap between the valence and the conduction band in the bilayer graphene. Calculating the 
occupation of the impurity level and the susceptibility in the bilayer graphene we show that 
the magnetic moment decreases with increasing the inter-layer coupling. 



2 



II. BILAYER GRAPHENE 



Fig{T] shows the lattice structure of the bilayer graphene with the adatom. The inter-layer 
stacking is assumed to be the Bernal order where the top layer has its A sublattice atop 
the sublattice B of the bottom layer. The bias voltage V is applied across the layers. The 
system Hamiltonian 

H = Htb + H, + Hf (1) 
contains the graphene bilayer term Htb, which in a tight-binding approximation reads 

2 

Htb = T.Hi + Ht + Hv, (2) 

1=1 

with 

Hi = -tY. [almbUHj] + 6L(R,)a,.(R,)], (3) 
Ht = -tpY.[aURi)b2a{Ki) + 6L(Ri)ai,(R,)], (4) 

Hv = ^ ^[aL(Ri)ai,(Ri) + 6L(R,)6i,(R,) - aL(Ri)a2<.(R.) - 6L(Ri)-&2.(R^)], (5) 

The operator a/o-(Ri) {bicr(R.i)) annihilates a state with a spin a at the position Rj on the 
sublattice A(B) of the / plane, t is the nearest neighbour in-plane hopping energy, tp is the 
inter-layer hopping energy. For the hybridization with the localized impurity states we write 

H. = VfJ2[f%A0) + bl{0)U (6) 

cr 

where faifl) is the annihilation (creation) operator of a state with a spin a at the impurity, 
and Vf is the hybridization strength. In the momentum space we have 



Hi = -tY^imalM^ + r(k)6La,k.], (7) 

ko- 

HT = tpY^ [aik<x^2ka + 4ka«lka] , (8) 
k,cr 

= IZ[aIka«iko + &Iko^ik<T - aLa«2ka - 4k<7^2ka], (9) 

kcr 

^^ = -^E(/i^ik. + 4./.), (10) 

V A 
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where 0(q) = -t E e'^^-^' with 5^ = f (1, v^, 0), 62 = f (1, -v^, 0), ^3 = a{l, 0, 0) (here a is 

i=l 

the lattice spacing), and N is the number of sites on sublattice B of plane 1. Diagonalizing 
Htb we find the spectrum 



E^^ik) = ±^el + j + ^±^^ + {tl + V-)el (11) 

where efc — ±|(;6(k)| is linearizable around the K points of the Brillouin zone by = ±7;^|k| 
where vp = 3ta/2 is the Fermi velocity. The impurity is described by Hamiltonian Hf with 

Hf^J2^ofja + Un^n^, (12) 

where — fife is the occupation number operator, £0 is the single electron energy at the 
impurity. The Coulomb interaction is included as a finite Anderson term U . For simplicity, 
we adopt a mean field approximation to the electronic correlations at the impurity, Un-irii ~ 
U Y^in-w) flfa — U{ni^){ni). The impurity Hamiltonian is rewritten as Hf = Y^^aflfa with 

a a 

£^(T = ^0 + U{n-^). To investigate the localized magnetic states, we calculate the occupation 
number of the electrons of a given spin a at the impurity. At low temperatures all the 
states below the Fermi level /x are completely occupied and the occupation of the impurity 
is determined by 

M = r dujp%u), (13) 

J —00 

where Pf{uj) is DOS at the impurity level. We infer it from the retarded Green's function 

G^r'it - 1') = -i^it - mfM fim- (14) 



Gr (^) = , (15) 



By the standard equation of motion, we can derive 

i 

cu — Ea — Ey-(a;) + ir]^ 
where 

•^^^ ^ t 4lk|^-2(-2 + l^)4|k|2 + (c.^-^)2-t2(^2_Zi)- l-LDJ 

Introducing a high-energy cutoff D of the graphene bandwidth, we obtain for uj > 

+^#{ ";;'^:^g!:rj^^ [sgn(t^)^(0 < < i^^) - sgn(^)^(0 < X. < D^)] 



-^^[sgn(f^)^(0 < < + sgn(^)^(0 < < D^)]], 



(17) 



where 9{x) is the step function, and 
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xi,2 = {oo' + ^)± 2 V4(^' + t>' - tiy- (18) 



For u^<',M 



4 tj+V^ ■ 



^rf ^ _ r u,-V/2 1 I g4_2(a;2+y2/4)j)2^(^2_y2/2)2_t2(^2„y2/4) ^ 

^/l^'' - D2[ 2 ^"^1 (c^2_y2/2)2_t2(^2_y2/4) | 

-^^^^(^^^^H±iV^(arctan . .^r-^-^.^/^ . + arctan ^^±^ 



^i2V2/4-(V2+t2)a;2 ^ y/tlV^/4-{V^+t2)uj2 ^t2y2/4_(v2+j2)^2^J- 

(19) 

The summation over k in Eq. (fT6l) is accurate for u D hj ensuring the conservation 
of the total number of states in the Brillouin zone according to the Debye's prescription. 
Substituting Sj(^) Eg. (113]) . the retarded Green's function G'j'^'^(co') can be obtained. 
Note, the determination of (n^) in Eq. flTSl) entails a self-consistent calculation of DOS at 
the impurity level via the relation p'f{uj) = — ^lmG'j'^'^(co'). When tp = V = 0, our present 
results reduce to those of Ref. [23 1. 



III. TRILAYER GRAPHENE 

The Hamiltonian for trilayer graphene contains a coupling the B atom of the second layer 
to the A atom of the third layer according to the conventional Bernal-type stacking order. 
Similar to the bilayer graphene case we find for the impurity Green's function 

to ~ + ir]'' 

where 



Gr (^) = . . , (20) 



yr _ Ai4|k|4+i?i4|k|2+Ci , . 

^ ~ N ^ l,^|k|6+B24|k|4+C2f^|k|2+D2 V^-^^ 

With = uj-^, B, = tluj-{uj-^)[ij^ + iij + ^)% C, = LU^{LU~^)iij + ^y-2tluj\u; + ^), 
B, = -3u^-^, C2 = -2ty + 3u^ + ^, D2 = -uj\uj^-y^Y + 2ty{uj^-y^). Performing 
the summation over k in Eq.(l2ll) as Eq.([l6l) we find for A = (25|-9fi2C2+27L'2)^+4(-5| + 
3^2)3 > the resuh 



1)2 1 [ 1 (x2—xi)(xs,—x\) J I X\ 



+ ^M£2±£ll±£^(arctan -^"^+"^^/^ + arctan , ^"^+"^^/^ ) 

y X2X3-(x2+x-i)'^ /A \/x2Xz~{x2+X3)'^ /i y/x2X-i-{x2+xzy /i 

, A]^{x2+X3)xi+BiXi^-AiX2Xs,+Ci r -1 D'^-{x2+X3)+X2Xz (22) 

{X2—Xl){x3—X\) I- 2 X2X3 

+ ^£i±£ilte^ (arctan ,^'-("^+"^)/^ + arctan , ("^+"^)/^ )]} 

2^J X2X-i-{x2+xz)'^ /A y X2Xz-{x2+xzY /A yj X2X-i-(x2+xz)'^ / 4: 

-zsgn(t^)^(0 <x,< D-Y-^[A, + -^-^--^t"S:;=ff^"^^^ ]' 
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where 

a^i = -f + ^l{-2Bl + 952^2 - 27/^2 + y'(25f - 952^2 + 27^2)' + 4(-5| + 3^2)3}^ 
+ -^{-2Bl + 952C2 - 27^2 - ^{2BI - 9B2C2 + 27^2)2 + 4(-5i + 302^}^ 

(23) 

X2,3 = -f + ^^i^^{-2Bl + 952^2 - 27^2 ± ^(251 - 952C2 + 27^2)2 + 4(-fi| + 3^2)3}^ 



I 2^' 2 1 



^{-253 + 9S2C2 - 27/^2 T ^(2fif - 952^2 + 27^2)2 + 4(-52 + 3^2)3}^- 

(24) 



For A < 0, 



_ _YLi\ A I [^l(3:^2+3:^3)+-Bi]3:^i-AiX2X3+Ci -1 1 I D^-xi I 
- £,2t[^l-h (a,.2-a;i)(a;3-X'i) x'l I 

r Ai(x2+3:3)+-Bi _|_ [Ai(x2+X3)+Bi]xi-AiX2Xs+Ci -\ i D^-X2 I 

I t Ai{x2+X3)+Bi I [Ai(a-2+X3)+Bi]a:i-Aia:2X3+Ci _ [Ai{x2+X3)+Bi]xi-AiX2X3+Ci i i -D^-X3 i -j 
I- X3-X2 (a;2-a;'i){x3-X2) (x2-a;i)(x3-zi) J I x'3 ' J 

_Sgn(fa)^(0 < X2 < D2)rAi(x2±£3)±Bi ^ [Ai(x2+x.3)+Bi]xi-AiX2X3+C7i 1 
^ ^ duj ' ^ ^ 'L a;3_2;2 \X2—x\)[x^—X2) J 

_)_Sp;n('^^3"\^('g < 3; < ^ A\(x2+xz)+B\ _|_ [Ai(a;2+a;3)+_Bi]3;i-AiX2X3+Ci _ [Ai(a;2+X3)+Bi]zi-Ai3;2X3+Ci i -j 

\ /V o /L x-i—X2 (x2—X\)(x-i,—X2) (x2—X\){x-3,—X\) JJ ' 

(25) 

where 

= + V^|-3C2 ^Qg(™T)^ = -f + ^V^ME^^ C0S( 2-+™^ ), (26) 



^ _ B2 I V-BI-3C2 ^^„/ 47r+arccosT N y _ 2(i3|-3C2)g2-3(B2C2-9D2) f,^^-. 
^3 3 3 ^ 2(B|-3C2)^ ■ ^ 

Substituting Eqs. fl22|) and (!25l) in Eq. (!20|) . we can derive self-consistently the occupation on 
the impurity for case of a trilayer graphene. 



IV. NUMERICAL ANALYSIS 



From the occupation of the two spin channel on the impurity we conclude on the formation 
of localized magnetic moment whenever 7^ nj^. For a detailed study conventionally, one 
introduces the dimensionless parameters 

X = DT/U and y = (i^ - eo)/U with F = nV^/D\ (28) 

The transition curves from the magnetic to the non-magnetic behavior as a function of the 
parameters x and y for the different hybridization and inter-layer coupling in the bilayer 
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graphene are shown in Figj2l For tp = V = 0, our results reduce to those of Ref.[23j: The 
magnetic boundary exhibits an asymmetry around y = 0.5, and can even cross when hne 
y = I. The magnetic region shrinks in the x direction with the hybridization Vf is increased; 
for y close to 1 (cf. eq. fl25]) ). the boundary line for magnetic transition shifts away from the y 
axis due to the increased influence of graphene on the impurity magnetization with enhanced 
hybridization. When the inter-layer coupling tp is taken into account (see Figj2](b)), the size 
of the magnetic region diminishes rapidly, and for a large enough tp, the magnetic boundary 
shrinks above the line y = 0. However, the magnetic boundary does not turn symmetric 
around y = 0.5, and the above magnetic boundary line crosses the line y = 1. The origin 
of this phenomena lies in the peculiar nature of the quasiparticles in the bilayer graphene; 
they exhibits features akin both to Dirac and to conventional fermions. The contribution 
of conventional fermions originates from the interlayer coupling that supports a metallic 
bilayer graphene and results in effects as for a conventional metallic host on the magnetic 
properties of the impurity. For large interlayer coupling we observe therefore magnetic 
boundaries similar an impurity in an ordinary metal. (FigOI) shows for a bilayer graphene 
the boundary between magnetic and non-magnetic impurity states as a function of the 
parameters x and y (eq j28|) for different impurity energy levels ^o- For V = the size of 
the magnetic region grows as Eq approaches the energy of the Dirac point. This behavior is 



reminiscent of the single layer of graphene 



23|, and originates from the suppression of the 



DOS around the impurity energy level. In contrast, for a nonzero gate bias, when Eq is 
close to the Dirac point from the positive energy side, the size of the region first increases 
to the maximum, then decreases with decreasing Eq, as shown in Figl3](b). The explanation 
for this phenomenon is as follows: the gate bias voltage gives rise to a finite electronic gap 
between the conduction and the valence band, and induces a large local DOS close to the 



gap edge s [24l| . In particular, the DOS may extend into the gap due to the influence of the 
impurity {25]. In this situation, the coupling between the bath and the impurity is enhanced 
inside the gap as compared with the zero bias case, leading thus to the non-monotonic 
dependence of the size of the region with Eq- 

FigJH shows the magnetic transition curve as a function of the parameters x and y (eq j28|) 
for different Eq in the trilayer graphene. For = 0, phenomena such as the asymmetry 
around the line y = 0.5 and the crossing of the line y = 1 in the magnetic boundary suggest 
the existence of Dirac fermions in the trilayer graphene. As eq approaches the energy of the 
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Dirac point, the magnetization region of the impurity grows due to the two almost-hnear 



touched bands reminiscent of the bands in monolayer graphene|26l|. It is interesting to note 
that for nonzero gate bias, the impurity magnetization region increases monotonously when 
Eq is close to the Dirac point, which is clearly different from that in the bilayer graphene. This 
behavior stems from the fact that the gate bias can not destroy the particle-hole degeneracy 
in the trilayer graphene |26|. 

To investigate the localized magnetic moment of the impurity in the magnetic region and 
the magnetic transition we calculate the magnetic susceptibility. The energy of the impurity 
spin states in a magnetic field B is = £o — crfJ^sB + Uria. The magnetic susceptibility of 
the impurity derives from 

de— dsa- 

FigjS] shows the occupation of the impurity spin level and the magnetic susceptibility as a 
function of y for the different inter-layer coupling in a bilayer graphene. The occupation (no-) 
versus y is a bubble that corresponds to the impurity magnetization. The corresponding 
susceptibility exhibits two peaks at the magnetization edge indicating the strength of the 
magnetic transition. For tp = 0, a strong magnetic moment of ~ 0.7 forms in almost 
the whole magnetic region. With increasing the inter-layer coupling tp, the magnetic bubble 
region diminishes signalling the decrease of the magnetic moment of the impurity, and the 
magnetic transition becomes very sharp. There is no localized magnetic moment in the case 
of a sufficiently strong inter-layer coupling. In this case, the magnetic boundary shrinks 
below the line a; = 6 in the x direction(see Figl2](b)). Figl6] shows the occupation of the 
impurity level and the magnetic susceptibility as a function of y for the different impurity 
energy level Eq in the bilayer graphene. The corresponding magnetic boundaries are defined 
in Figl3] (a) and (b) respectively. For V = 0, the magnetic bubble shifts towards the (n^) 
axis, and decreases with increasing Eq. When Eq becomes large enough, the bubble vanishes, 
meaning that the impurity loses magnetism in this situation. For large Eq the magnetic 
transition becomes very sharp. Inspecting FigEl(c) and (d) we find when the gate bias V is 
applied, the magnetic bubble shows a non-monotonic dependence on Eq, while the magnetic 
transition becomes very sharp with increasing Eo- Since the magnetic boundary line shrinks 
in the left hand side of the line x = 4.2 at Eq/D = 0.082(see Fig|3](b)), the impurity remains 
non- magnetic for any y, i.e. n| = n^, as shown in Figj6] (c). 
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V. CONCLUSIONS 



Summarizing, we studied the localized magnetic states of an impurity in biased bilayer and 
trilayer graphene. We find that the size of the magnetic region decreases rapidly compared 
with that in monolayer graphene, the impurity can magnetize even when the energy of the 
doubly occupied state is below the Fermi level, and the impurity magnetization region has 
a different shape. We can trace this behaviour back to the special nature of quasiparticles. 
When a gate bias is applied, the dependence of the impurity magnetic states on the impurity 
energy for the bilayer graphene shows a behavior different from that for a trilayer graphene 
due to the opening of a gap between the valence and the conduction band in the bilayer 
graphene. Correspondingly, the magnetic moment of the impurity versus the impurity energy 
in the bilayer graphene is affected strongly by the band gap induced by the gate bias. 
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FIG. 2: Regions of the magnetic and the nonmagnetic phase for the bilayer graphene. The boundary 
line gives y as a function of x (defined in eq. (|28p ) at tp/D = for the different Vf/D (a), and at 
Vf/D = 0.14 for different tp/D (b). The other parameters are Eq/D = 0.029 and V/D = 0. 
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FIG. 3: Regions of magnetic and non-magnetic phase for the bilayer graphene. The boundary Une 
gives y as a function of x for the different eo/D aX V/D = (a) and dX V/D = 0.05 (b), where 
Vf/D = 0.3 and tp/D = 0.05. 
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FIG. 4: Regions of magnetic and non-magnetic phase for the trilayer graphene. The boundary hne 
gives y as a function of x for the different £q/D at V/D = (a) and at V/D = 0.05 (b), where 
Vf/D = 0.2 and tp/D = 0.05. 
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1.0 




FIG. 5: The occupation of the impurity spin level and the magnetic susceptibility in the bilayer 
graphene for the different tp/D at x = 6. The other parameters are those of Figl2]^b). 
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FIG. 6: The occupation of the impurity spin level and the magnetic susceptibility in the bilayer 
graphene for different Eq/D at V/D = and x = 3.2 (a)-(b), and at V/D = 0.05 and x = 4.2 
(c)-(d). The other parameters are the same as in Figl3] 
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